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ABSTRACT 

We estimate the distribution of intrinsic shapes of UZC-SSRS2 groups of galaxies from 
the distribution of their apparent shapes. We measure the projected group axial ratio 
using the moments of their discrete galaxy distribution. Then using the non-parametric 
kernel method to estimate the smooth apparent axial ratio distribution we numerically 
invert a set of integral equations to recover the corresponding intrinsic distribution 
under the assumption that groups are either oblate or prolate spheroids. We find that 
the prolate spheroidal model fits very well the UZC-SSRS2 group distribution with a 
true mean axial ratio (/3) ~ 0.3 and up ~ 0.15. This shows that groups of galaxies 
are significantly more elongated, both on the plane of the sky and in 3 dimensions, 
than clusters of galaxies. The poorest groups that we consider, those with 4 members, 
are even more elongated than the overall population with 85% of the groups having 
I3< 0.4. 
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1 INTRODUCTION 



In Cold Dark Matter models, structure formation evolves 
in a hierarchical fashion with aggregation of smaller mass 
units along large-scale anisotropic structures. Since virial- 
ization processes tend to sphericalize initial anisotropic dis- 
tributions of matter, the shape of cosmic structures is related 
to their formation processes and evolutionary stage and thus 
it is extremely important to unambiguously determine their 
intrinsic shapes. 

Apart from disk galaxies all cosmic structures on larger 
scales appear to be dominated by prolate like shapes. This 
has been shown to be the case for clusters of galaxies (cf. 
Carter & Metcalfe 1980; Plionis, Barrow & Frenk 1991; 
Cooray 1999; Basilakos, Plionis & Maddox 2000) as well 
as for superclusters which show a predominance of filamen- 
tary like shapes both observationally, theoretically and in 
Cosmological N-body simulations (cf Zeldovich, Einasto & 
Shandarin 1982; Shandarin & Zeldovich 1983; Broadhurst et 
al. 1990; de Lapparent, Geller & Huchra 1991; Plionis, Jing 
& Valdarnini 1992; Jaaniste et al 1998; Sathyaprakash et al. 
1998; Valdarnini, Ghizzardi & Bonometto 1999; Basilakos, 
Plionis & Rowan- Robinson 2001). In the case of Hickson or 
Shakhbazian compact groups it has been shown that they 
are even flatter than clusters and most probably prolate-like 
configurations with typical true axial ratios ~ 0.3 (Varda- 
nian & Melik-Alaverdian 1978; Hickson et al. 1984; Malykh 
& Orlov 1986; Oleak et al. 1995). 

It is obvious that the intrinsic shape of cosmic struc- 
tures can be lost when projected on the plane of the sky 



and therefore it is important to deal with such and other 
systematic effects that can hide the true shape of cosmic 
structures. Different studies have attempted to recover the 
distribution of intrinsic shapes from the corresponding ap- 
parent distribution using inversion techniques based on the 
assumption that their orientations are random. 

The plan of the paper is the following: In Section 2 we 
describe the group sample that we use. We attempt to iden- 
tify the extent of projection contamination and we describe 
the projected shape determination method. In Section 3 we 
invert the projected axial ratio distribution and recover the 
corresponding intrinsic one. The discussion and our conclu- 
sions are presented in Section 4. 



2 PROJECTED POOR GROUP SHAPES 

We use the recent UZC-SSRS2 group catalogue (Ramella et 
al 2002) which is based on the Updated Zwicky Catalogue 
(UZC; Falco et al. 1999) and the Southern Sky Redshift Sur- 
vey (SSRS2; da Costa et al 1998), to measure the projected 
group shape distribution and hence attempt to estimate 
their intrinsic shape. The catalogue has a sky coverage of 
4.69 sr, it is limited to ms ~ 15.5 and contains 1168 groups 
of galaxies exceeding a number density contrast threshold 
of 5p/p = 80. The group catalogue was constructed using 
the well-known friends of friends algorithm with a linking 
parameter that scales with increasing redshift, in order to 
take into account the galaxy selection function. 
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Figure 1. The group number density in equal volume shells. Note 
that we plot the number density only for cz > 1800 km s^^. At 
smaller distances there is a local excess in detected groups. The 
dashed line shows a crude fit to the data 
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Figure 2. The dependence of the group velocity dispersion (left 
panel) and maximum intergalaxy separation in Mpc (right 
panel) on the group redshift. The continuous lines represent the 
dependence of the line-of-sight velocity (Vl; left panel) and the 
projected separation (D^ ; right panel) galaxy pair linking param- 
eters on the redshift. 



In order to have a representative sample of the true un- 
derlying group population, we have chosen to study those 
groups the number density of which, within some limiting 
redshift, is relatively constant. We have derived the group 
number density as a function of redshift in equal volume 
shells (see Fig. 1) and we have found that it is roughly con- 
stant out to cz ~ 5500 km s~^, which we choose as our limit. 
Within this velocity limit we have a total of 245 groups with 
4 or more members. 

In Fig. 2 we present the group's velocity dispersion as 
well as their maximum intergalaxy separation as a function 
of redshift. A strong ^-dependence is evident, although such 
dependence is rather weak within ~ 5500 km s~^. What is 
the cause of this redshift trend ? 

We note that in order for the Ramella et al. (2002) al- 
gorithm to identify groups having the same limiting density 
contrast at the different distances, it was necessary to take 
into account the drop with redshift of the galaxy space den- 
sity which is due to the magnitude limit of the parent galaxy 
catalogue. The way this is accomplished is by increasing with 
redshift the linking parameters, used to identify the group 
members. The Ramella et al (2002) algorithm, like most 



others, uses two linking parameters; one that links pairs of 
galaxies below some projected separation, ie., D12 < Dl, 
and one that links galaxies with line-of-sight velocity differ- 
ences below so threshold, ie., Vi2(= |Vi — V2I) < Vl- There- 
fore in order to keep the limiting density enhancement of the 
detected groups constant it is necessary to scale the linking 
parameters by a distant dependent quantity that compen- 
sates, as discussed before, for the drop of the selection func- 
tion of the parent galaxy population. For example, Ramella 
et al (2002) use the following scaling: 



Dl = DoR and Vl = VoR , 



where 



R 



<i>{L)dL/ 



<i>{L)dL 



1/3 



with $(L) the galaxy luminosity function, Lmin the faintest 
luminosity at which galaxies with the limiting magnitude 
of the catalogue can be visible at the fiducial velocity used 
(1000 km/sec) and Lmin{r) the faintest luminosity at which 
galaxies with the limiting magnitude of the catalogue can 
be visible at the distance r. In Figure 2 we also plot by 
continuous lines the values of Vl (multiplied by a factor 
1/5) and Dl (multiplied by 1.5) in the corresponding plots, 
which evidently follow the redshift trend observed in the 
data. 

The apparently good correlation between the increase 
with redshift of the size and velocity dispersion of groups 
and the corresponding increase of the group linking param- 
eters suggests that the bias is introduced by the algorithm 
itself, especially by the way it deals with the drop of the 
galaxy redshift selection function. The fact that groups of 
galaxies are clustered (eg. Zandivarez, Merchan & Padilla 
2003 and references therein) implies that the increase with 
redshift of the group linking parameters may increasingly 
include, as part of the identified groups, galaxies belonging 
to neighbouring structures; a fact which will affect both the 
size and velocity dispersion of the resulting groups. This sys- 
tematic bias could be quantified with the use of simulations 
but it is out of the scope of the present work. 

We conclude that the probability that the groups identi- 
fied constitute a homogeneous and biased free sample as well 
as that they constitute real dynamical entities, decreases 
with redshift. 



2.1 Shape determination methods 

We have used two methods to determine the group shapes. 
The first one is according to Rood (1979), by which the axial 
ratio, q = (b/a), is such that a is the angular distance be- 
tween the most widely separated galaxies in the group, and b 
is the the sum of the angular distances 61 and 62 of the most 
distant galaxies on either side of the line a joining the most 
separated galaxies. The second one is the moments of inertia 
method (cf. Basilakos et al 2000 and references therein). In 
Fig. 3 we compare the results of the two methods. It is ev- 
ident that they both provide equivalent axial ratios, except 
in a few cases where the position of the largest pair sepa- 
ration, used to determine a in the first method, is strongly 
asymmetrical with respect to the rest of the galaxies of the 
group (the outliers in the figure). 
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Figure 3. Comparison of tiie group axial ratio determined using 
the moments of inertia and Rood's (1979) methods. 

In what follows we will use the results of the more ro- 
bust moments of inertia method. Due the the slightly more 
complicated nature of this method, we describe it in more 
details below. Firstly, the galaxy equatorial positions are 
transformed into an equal area coordinate system, centered 
on the group center of mass, using: x — (ug — agr) x cos(5gr) 
and y = Sg~ Sgr, where subscripts g and gr refer to galaxies 
and the groups, respectively. We then evaluate the moments: 

hi = - x^) 

i 

I22 = y^^Wijri -yf) 

i 

I12 = I21 = —^'^■WiXiyi (1) 

i 

with Wi the statistical weight of each point (in our case Wi = 
1) and Ti the distance of the j"* galaxy from the group center 
of mass. Note that because the inertia tensor is symmetric 
we have /12 — hi- Diagonalizing the inertia tensor 

det(Jij - = (M2 is 2 X 2 unit matrix.) (2) 

we obtain the eigenvalues Ai, A2, from which we define the 
principal axial ratio of the configuration under study by: q = 
A2/A1 = b/a, with Ai > A2. The corresponding eigenvectors 
provide the direction of the principal axes. 



2.2 Projection effects 

Random projections of field galaxies could appear as groups 
of galaxies, which however would bare no relation to dy- 
namical entities. We have performed Monte-Carlo simula- 
tions {Nsim = 10000) and found that sets of 4 points uni- 
formly distributed within a sphere and projected on a plane 
give a nearly Gaussian g-distribution with (q) ~ 0.6 and 
(Tq ~ 0.18, while for sets of 10 points the corresponding val- 
ues are (q) ~ 0.72 and aq ~ 0.13 respectively, significantly 
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Figure 4. Comparison of the axial ratio distribution of the groups 
with 4 members (continuous line) with the corresponding distri- 
bution of random Monte-Carlo groups (dashed line). 



larger than the observed case (see Table 1). As an exam- 
ple we compare in Fig. 4 the axial ratio distribution of our 
groups (having 4 galaxy members) with the corresponding 
Monte-Carlo random realizations. The difference is evident. 
However, this does not preclude the possibility of projec- 
tion effects altering the apparent shape and the dynamical 
parameters of real groups. 

Ramella, Pisani & Geller (1997) and Diaferio et al. 
(1999) based on either geometrical or N-body simulations 
and reproducing the group identification procedure have 
found that between 70% and 80% of the UZC-SSRS2 groups 
should be real dynamical entities, while the rest are expected 
to be the result of superpositions of field galaxies. Such a 
fraction of false groups, however, could be a serious prob- 
lem for studies of the group dynamical properties, and thus 
should be investigated in detail. 

Indeed, during our eye-balling process of verifying the 
group shape measurements, we encountered some cases with 
"strange" positional and velocity configurations, which ap- 
peared indeed to be the results of projection effects. This 
prompt us to check the most probable candidates of be- 
ing fake groups, those with a relatively high-velocity disper- 
sion. One such example is the group UZC269, containing 
four galaxies, three of which (forming a straight line) have 
{v) ~ 4450 km s~^ and velocity dispersion — 200 km 
s^^, while the fourth galaxy has a velocity 4ct„ lower than 
(v). Another example is UZC156 with 6 members and an 
overall (j„ — 647. However, the total projected size of the 
group is ~ 2.0 Mpc while it can be separated in both veloc- 
ity space and on the plane of the sky, in two triplets having 
(v) = 5200±290 km s"^ and (v) = 4200±195 km s~\ Other 
fake cases appeared to contain two independent groups or 
pairs of galaxies, with large velocity and/or spatial separa- 
tions. 

Of course the probability of a group being false is in- 
versely proportional to the group galaxy membership, rim- 
Random projections will affect significantly more the appar- 
ent characteristics (dynamical and morphological) of small 
groups rather than large ones and for this reason we have 
decided to exclude from our study groups with rim ~ 3, ex- 
actly due to their a priori high probability of being chance 
projections and not real physical systems (see also Focardi 
& Kelm 2002). 
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Table 1. The median axial ratio, q, and the median value of the 
maximum intergalaxy separation, a (in Mpc), of the differ- 
ent membership UZGC groups with cz < 5500 km (after 
exclusion of candidate false groups). 
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Furthermore, we have devised an objective algorithm 
in an attempt to single out candidate false groups having 
rim = 4 and 5 (which are those affected the most). We will 
present this algorithm in an accompanying paper, which 
deals with the dynamics of groups (Tovmassian & Plionis 
2003 in preparation). Here we only present the basic as- 
sumption on which the algorithm is based, which is that the 
probability of a group being false increases with increasing 
deviation of its member velocities from a Gaussian distribu- 
tion having mean and variance the observed group values. 
We find in total 20 and 15 groups with n-m = 4 or 5, re- 
spectively that are probably the result of superpositions of 
galax;ies. Note that these groups constitute a relatively small 
fraction (~ 14%) of the sample of 245 groups, less than the 
expected number of false groups, according to Ramella et 
al. (2002). 

We have also tested our results (presented in section 3) 
by including these groups and found that indeed they do 
create problems in the inversion from the projected to the 
intrinsic 3D axial ratio distribution (they cause the inverted 
distribution for the prolate case to have unphysical negative 
values; see section 3). 

Since the probability of a group being affected by pro- 
jections is inversely proportional to Um, we will present re- 
sults also separately for groups with Um = 4 and 4 < nm < 
10. 



2.3 Mean group shape peirameters 

The summary of the main structure parameters of the dif- 
ferent membership samples of groups is presented in Table 
1. The first and second columns give the group membership, 
rim, and the number N of such groups, respectively, while 
the third and fourth columns show the median q and a (ma- 
jor axis) values, together with their 68% and 32% quantile 
values. It is evident that the considered UZC-SSRS2 groups 
are very elongated, significantly more than what expected 
from random projections of field galaxies, giving support 
to them being real dynamical entities. Furthermore, among 
these groups there are many so called chain-like groups, with 
q smaller than 0.20-0.30. The large number of chain-like 
groups, which also determines the relatively small median 
values of q, is similar to that of compact groups, the pre- 
dominance of which was first mentioned by Arp (1973). 

Note that the median q of all poor groups (4 < rim < 
10) is 0.33^g Qg. A certain correlation is apparent with the 
median b/a and a increasing with nm- The increase of the 
major axis with nm could have been due to the increase with 
redshift of the group linking parameter (see Ramella et al. 



2002), a fact that induces the systematic increase with z of 
the group size (see Fig. 1). However, we have verified that 
the mean redshift (see Table 1) of each group subsample is 
constant and thus the above systematic effect is apparently 
not the cause of the observed increase of a with z. 

The increase of the group sphericity could be explained 
as an indication of a higher degree of virialization, which 
is expected to be more rapid in systems containing more 
galaxies (mass). However, the increase of their major axis, 
if proven to be true, is somewhat more perplexing. 

2.4 The projected axial ratio distribution 

We now proceed in describing our approach of fitting the 
observed discrete distribution of axial ratios with the so- 
called kernel estimators. Here we review the basic steps of 
the Kernel method, following the notation of Ryden (1996) 
but for further extensive reviews see Silverman (1986), Scott 
(1992), Vio et al. (1994) and Tremblay & Merritt (1995). 

Given the sample of axial ratios qi,q2, ■■■■,qN for N 
groups, the kernel estimate of the frequency distribution is 
defined as: 

N 
i 

where K{t) is the kernel function, defined so that 

/ + 00 
K{t)dt = 1 , (4) 
- oo 

and h is the "kernel width" which determines the balance 
between smoothing and noise in the estimated distribution. 
In general the value of h is chosen so that the expected value 
of the integrated mean square error between the true, f{q), 
and estimated, /(g), distributions, [/if (a;) — f{x)\ dx, 
is minimised (cf. Vio et al. 1994; Tremblay & Merritt 1995). 
As it has been shown in different studies the choice of a 
kernel function, K{t), among quadratic, quartic and Gaus- 
sian forms (cf. Tremblay & Merritt 1995), provide fits that 
differ trivially only in their asymptotic efficiencies. We have 
chosen a Gaussian kernel: 

m = ^e-*'/' . (5) 
V 27r 

In order to obtain physically physically acceptable re- 
sults with /(g) = for g < and g > 1, we apply reflective 
boundary conditions (cf. Silverman 1986; Ryden 1996), re- 
placing the Gaussian kernel with: 

K„„...) = «(i^).A-(l±*)+.-(^) 

which also ensures the correct normalization, f[q)dq= 1. 

In Figure 5 we present the projected axial ratio dis- 
tributions for the different membership groups (circles), as 
indicated in the different panels, with their Poisson la error 
bars, while the solid lines shows the kernel estimate / for 
the appropriate width, h. 

3 TRUE GROUP SHAPES 

In order to find the intrinsic axial ratio distribution assum- 
ing that groups are either oblate or prolate spheroids, we 



© 0000 RAS, MNRAS 000, 000-000 



Shapes of Poor Groups 5 




Figure 5. The apparent axial ratio distributions for different 
group membership. The sohd hne is the smooth fit from the non- 
parametric kernel estimator. 
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Figure 6. The distribution of the intrinsic poor group axial ra- 
tios (continuous line) and its la range (broken lines) assuming 
that they are either prolate (left panel) or oblate (right panel) 
spheroids. 



use an inversion method, described below. Although there is 
no physical justification for the restriction to oblate or pro- 
late spheroids, it greatly simplifies the inversion problem. 
Furthermore, if groups are a mixture of the two spheroidal 
populations or they have triaxial configurations then there 
is no unique inversion (Plionis, Barrow & Frenk 1991). 

The relation between the apparent and intrinsic axial 
ratios, is described by a set of integral equations first investi- 
gated by Hubble (1926). These are based on the assumptions 
that the orientations are random with respect to the line of 
sight, and that the intrinsic shapes can be approximated by 
either oblate or prolate spheroids. Writing the intrinsic ax- 
ial ratios as /3 and the estimated distribution function as 
No{P) for oblate spheroids, and Np{f3) for prolate spheroids 
then the corresponding distribution of apparent axial ratios 
is given for the oblate case by: 



No{l3)d(3 



to (l-g2)l/2(g2_ ^2)1/2 

and for the prolate case by: 



/327Vp(/3)d/3 



g2 (I_q2)l/2(g2_ ^2)1/2 



(7) 



(8) 



Inverting equations (eq.7) and (eq.8) gives us the distribu- 
tion of real axial ratios as a function of the measured distri- 
bution: 

2/3(1^/32)1/2 1-0 



and 

NAP) = 



2(1-/3^)1/2 ^^^^^^^ 



n/3 7q dg''' ''(/ja -q2)i/2 



(9) 



(10) 



with /(O) = 0. In order for iVp(/3) and No{l3) to be phys- 
ically meaningful they should be positive for all /3's. Fol- 
lowing Ryden (1996), we numerically integrate eq.(9) and 
eq.(lO) allowing Np{P) and No{l3) to take any value. If the 
inverted distribution of axial ratios has significantly negative 
values, a fact which is unphysical, then this can be viewed 
as a strong indication that the particular spheroidal model 
is unacceptable. 

In Figure 6 we present the intrinsic group axial ratio 
distributions. The oblate model (right panel) is completely 
unacceptable since it produces only negative values of the 
inverted intrinsic axial ratio distribution. The UZC-SSRS2 
groups shape is represented well only by that of prolate 
spheroids which is in agreement with previous studies of al- 
beit smaller group samples (cf. Oleak et al 1995). It is very 
interesting the fact that there are almost no groups with 
true axial ratio, /3^ 0.6, ie., there are no roundish groups 
while most of them are extremely elongated. The most elon- 
gated groups are the poorer ones (those with rim = 4) with 
85% having /3< 0.4. We can crudely approximate the in- 
trinsic prolate axial ratio distribution of all the UZC-SSRS2 
groups of our sample (lower left panel of Fig. 6) by a Gaus- 
sian having {/3) = 0.29 and cr = 0.16 (see Fig. 7). 

How do our group shape results compare with those of 
clusters of galaxies ? In Figure 8 we present the distribution 
of intrinsic axial ratio, for the prolate case (which also is the 
best model for clusters; see Plionis et al 1991, Basilakos et al 
2000), for both the UZC-SSRS2 groups and the APM clus- 
ters that appear to have no substructure (which means that 
they are probably in virial equilibrium; see Basilakos et al 
2000). It is evident that groups are significantly more elon- 
gated than clusters and although there are quite a few spher- 
ical clusters, this is not the case for groups. This probably 
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Figure 7. The intrinsic axial ratio distribution (for the prolate 
case) of all the UZC-SSRS2 groups of the sample analysed (con- 
tinuous line) and the best Gaussian fit (broken line). 
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that generically poor groups of galaxies, compact or not, 
are extremely elongated prolate systems, much more than 
clusters or even elliptical galaxies. 

This result supports the view by which groups form by 
accretion of gala:xies along larger structures, like filaments, 
in which they could be embedded (cf. West 1994). Such 
an accretion process would happen preferentially along the 
group major axis, which should then be typically aligned 
with other nearby structures (groups or clusters). Such an 
effect is well documented for clusters (cf. Bingelli 1982; Plio- 
nis 1994, Plionis & Basilakos 2002 and references therein), 
but it has still to be determined for groups of galaxies. 

Finally, we would like to comment on the orientation ef- 
fects that the extremely elongated and intrinsically prolate 
nature of groups would create. When elongated groups are 
orientated roughly along the line of sight they will appear 
roughly spherical therefore having a higher velocity disper- 
sion and smaller sizes while elongated groups seen roughly 
orthogonal to the line of sight will appear to have smaller 
velocity dispersions. This would aflect relations like that be- 
tween the group's X-ray luminosity - velocity dispersion and 
thus if not taken into account, erroneous conclusions could 
be reached regarding their dynamical state (see Tovmassian, 
Yam & Tiersch 2002; Plionis & Tovmassian 2004). 
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Figure 8. Comparison of the UZC-SSRS2 group (continuous line) 
and the APM cluster (broken line; from Basilakos et al 2000) 
intrinsic axial ratio distributions. Note that these distributions 
are based on 212 and 405 objects respectively. 

implies that the cluster distribution is in a more advanced 
dynamical state, with violent relaxation having sphericalized 
the initial flattened distribution of galaxy members. 



4 DISCUSSION & CONCLUSIONS 

We have measured the projected axial ratio of all UZC- 
SSRS2 Groups of Galaxies within a volume-limited region 
{cz < 5500 km/sec) using the moments of the discrete galaxy 
distribution. 

Using the nonparametric kernel procedure we obtain a 
smooth estimate of the apparent UZC-SSRS2 group axial- 
ratio distribution. The projected axial ratio distribution of 
the whole group sample peaks at 5 ~ 0.33 with an extended 
tail towards apparently spherical groups. Assuming that the 
UZC-SSRS2 groups constitute a homogeneous population of 
either oblate or prolate spheroids, we numerically invert the 
apparent axial ratio distribution to obtain the corresponding 
intrinsic one. The only acceptable model is provided by that 
of prolate spheroids having an intrinsic distribution with 
(P) ~ 0.29 and = 0.16 (if modeled by a Gaussian). This 
results are in very good agreement with the analysis of 95 
Shakhbazian compact groups (Oleak et al 1995) and shows 
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